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ABSTRACT 

Modelling of adiabatic gravity wave propagation in the solar atmosphere showed that 
mode conversion to field guided acoustic waves or Alfven waves was possible in the 
presence of highly inclined magnetic fields. This work aims to extend the previous 
adiabatic study, exploring the consequences of radiative damping on the propagation 
and mode conversion of gravity waves in the solar atmosphere. We model gravity waves 
in a VAL-C atmosphere, subject to a uniform, and arbitrarily orientated magnetic 
field, using the Newton cooling approximation for radiatively damped propagation. 
The results indicate that the mode conversion pathways identified in the adiabatic 
study are maintained in the presence of damping. The wave energy fluxes are highly 
sensitive to the form of the height dependence of the radiative damping time. While 
simulations starting from 0.2 Mm result in modest flux attenuation compared to the 
adiabatic results, short damping times expected in the low photosphere effectively 
suppress gravity waves in simulations starting at the base of the photosphere. It is 
difficult to reconcile our results and observations of propagating gravity waves with 
significant energy flux at photospheric heights unless they are generated in situ, and 
even then, why they are observed to be propagating as low as 70 km where gravity 
waves should be radiatively overdamped. 
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1 INTRODUCTION 

Recent multi-height observations of low frequency oscilla- 
tions in the low solar chromosphere suggest that the energy 
flux carried by upward propagating gravity waves, with fre- 
quencies between 0.7 and 2.1 mHz, comfortably exceeds the 
co-spatial acoustic wave flux (Straus et al 2008). In our pre- 
vious paper (Newington & Cally 2010), (referred to hence- 
forth as paper I) , we explored the propagation, reflection and 
mode conversion of gravity waves in a VAL C related atmo- 
sphere, permeated by uniform, inclined magnetic fleld. The 
significant flnding of that study was that in regions of highly 
inclined magnetic field, gravity waves experience mode con- 
version to up-going (field-guided) acoustic or Alfven waves. 
While acoustic waves are likely to shock before reaching the 
upper chromosphere, Alfven waves can propagate to greater 
atmospheric heights, perhaps contributing to the observed 
coronal Alfvenic oscillations (De Pontieu et al. 2007; Tom- 
czyk et al. 2007). 

A limitation of the work presented in paper I was the 



* E-mail: Marie.Ncwington@inonash.edu 
t E-mail; Paul.Cally@monash.edu 



assumption of adiabatic wave propagation, which is known 
to be invalid in the photosphere and low chromosphere. A 
more realistic investigation of the propagation and mode 
conversion of gravity waves in the solar atmosphere would 
include the effects of radiative damping. 

Radiatively damped atmospheric gravity waves have 
been considered by a small number of authors previously, 
(see Souflrin 1966, Stix 1970, and the exhaustive study by 
Mihalas & Toomre 1982), but all of those studies have been 
purely hydrodynamic in nature, with no applied magnetic 
field, hence there is no possibility of mode conversion. Our 
principal objective in this work is to explore the conse- 
quences of radiative damping on the mode conversion of 
gravity waves, and this requires a magnetohydrodynamic 
(MHD) treatment. 

In this paper we extend our earlier MHD study, relaxing 
the assumption of adiabatic wave propagation by incorpo- 
rating radiative damping using the Newton cooling approx- 
imation. We recognise that this approximation is strictly 
valid only for optically thin perturbations of a homogeneous, 
infinite, and isothermal gas, and that none of these condi- 
tions are realised in the photosphere and chromosphere. The 
advantage of using this simple treatment is that the same 
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mathematical tools used in paper I can be employed with 
minor modifications. The insight gained from this simple 
model may direct our attention to cases that are interesting 
enough to warrant a more thorough treatment of the effects 
of radiation. 

The questions central to our investigation are: How are 
the mode conversion pathways affected by radiative damp- 
ing? Do we still get appreciable mode conversion to Alfven 
waves when damping is present? 

This paper is organised as follows: Section 2 presents the 
atmospheric model and the mathematical tools used in this 
investigation, namely the damped dispersion relation and 
numerical solution of the damped linear MHD wave equa- 
tions. Section 3 presents the results. Dispersion diagrams 
reveal insights into the mode conversion pathways in z-k^ 
phase space, and numerical integration of the wave equations 
quantify the acoustic and magnetic fluxes. The conclusions 
are stated in Section 4. 



2 MODEL AND EQUATIONS 

This section describes the atmospheric model, the coordi- 
nate system and the equations used to generate the results 
of this paper. As this work is an extension of the adiabatic 
investigation of paper I, the reader will be referred to that 
paper to avoid repetition of material, where appropriate. 



2.1 Atmospheric model and coordinate system 

We use the same atmospheric model as employed in paper I 
- the Schmitz & Fleck (2003) adaptation of the horizontally 
invariant VAL C model (Vernazza et al. 1981) up to a height 
1.6 Mm above the base of the photosphere. An isothermal 
top is appended above 1.6 Mm. No transition region is in- 
cluded. A uniform inclined magnetic field is imposed upon 
this atmosphere. 

The origin of the coordinate system is located at the 
base of the photosphere. The coordinate system is orientated 
such that the wavevector (fe) lies in the x-z plane, and the 
orientation of the magnetic field (B) is described in terms of 
the inclination from the vertical (9) and the azimuthal angle 
{(j)) (see Fig. 1). We distinguish between two dimensional 
(2D) and three dimensional (3D) configurations by the angle 
between the magnetic field and the vertical plane of wave 
propagation ((ji): 2D: </!> = 0; 3D: (j) ^ 0. 

Following paper I, our tools for investigation of radia- 
tively damped gravity wave propagation are the dispersion 
relation and numerical integration of the linearised MHD 
wave equations. Ray theory is not directly employed as it 
is difficult to institute and interpret in the presence of dis- 
sipation. The derivation of the appropriate forms of these 
equations including Newton cooling are described below. 




Figure 1. The coordinate system and geometry adopted for this 
analysis. B is the magnetic field vector and k is the wavevector. 



sumed a function of height z only. As r — > oo, the motions 
become adiabatic. The continuity equation and the momen- 
tum equation are unchanged in the Newton cooling approxi- 
mation, but the energy equation is modified by the inclusion 
of a term proportional to the temperature perturbation, as 
follows (see, for example, Cally 1984): 
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Here p is the density, p is the pressure, T is the temperature 
and c is the sound speed. The subscripts and 1 denote the 
background values and Eulerian perturbations respectively. 

Upon linearising equation (1), applying the WKBJ ap- 
proximation, and making use of the perfect gas law, the def- 
inition of the adiabatic sound speed c^ — ■ypo/po (where 7 
is the ratio of specific heats), and the magentohydyrostatic 
balance, the following equation for the pressure perturbation 
is obtained. 
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and g is the acceleration due to gravity. H the density scale 
height and ^ = (^, r), C,) is the displacement vector. In the 
adiabatic limit ^1 and ^2 reduce to 1. 



2.2 Mathematical tools incorporating Ne'wton 
cooling 

2.2.1 Energy equation and expression for the pressure 
perturbation 

In the Newton cooling approximation, the radiative damp- 
ing is parametrized by the radiative relaxation time (r), as- 



2.2.2 Height dependence of the radiative relaxation time r 

Although expressions for the radiative damping time have 
been provided by Spiegel (1957) and Stix (1970), (for con- 
tinuum and line emission, respectively), the actual values for 
the radiative relaxation time in the photosphere and chro- 
mosphere are only crudely known. 

In this paper it will suffice to adopt the simple linear 
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Figure 2. The behaviour of the damping ratio as a function of 
height for the VAL-C atmosphere and the Mihalas & Toomre 
(1982) prescription for the radiative damping times. 
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Figure 3. Height dependence of temperature in the model atmo- 
sphere used here. The VAL-C model is employed below 1.6 Mm ; 
an isothermal layer is appended above 1.6 Mm. 



radiative damping time introduced by Mihalas & Toomre 
(1982) to generate the results in that paper (curve 2): 

r (2) =50 + (2200/3)2, (5) 

where r is expressed in seconds and 2 in megameters. 

Very short damping times are known to suppress the 
propagation of gravity waves. When Newton cooling is in- 
cluded in the energy equation, the velocity of a vertically dis- 
placed fluid element may be described in terms of a damped 
harmonic oscillator, (see, Souffrin 1966 and Bray & Lough- 
head 1974). This allows identification of the damping ratio 
as \/{2'yNT), where A*' is the Brunt- Vaisala or buoyancy 
frequency. If the damping ratio is greater than one, the mo- 
tion is overdamped and the fluid parcel will return to the 
equilibrium position without oscillating. Applying this con- 
dition to their atmosphere, Mihalas & Toomre (1982) noted 
that the wave was overdamped below 0.1 Mm and so they 
started their simulations from this height. The atmosphere 
in this paper has a different height dependence of the Brunt- 
Vaisala frequency, and the location where the damping ratio 
is 1 is higher at about 0.2 Mm (see fig. 2). Hence, we expect 
gravity waves to be underdamped and propagating above 
0.2 Mm, and overdamped below this height. 

We also ran simulations using constant relaxation times 
from 200s to 1 ks to gauge the sensitivity of the wave prop- 
agation to the form of r, but (5) was used to generate all 
the figures in this paper. 

2.2.3 Dispersion relation 

The adiabatic 3D MHD dispersion relation described in 
Newington & Cally (2010) was derived using a Lagrangian 
formulation, which is not easily extended to nonadiabatic 
wave propagation. However, a dispersion relation for the 
nonadiabatic case can be readily constructed from the adia- 
batic relation, if the simplifying assumption of an isothermal 
atmosphere is adopted. This is not a bad approximation for 
the region in question (see figure 3). 

In an isothermal atmosphere, the density scale height 
is H = c^/jg, and so /ii = l. The equation for the Eulerian 
pressure perturbation, equation (2), then becomes, 



where, following Cally (1984), we define a new (complex) 
quantity ^ 

C = fl2C^ . (7) 

Because the other perturbation equations are un- 
changed in the Newton cooling approximation, this suggests 
that for an isothermal atmosphere, the equations describ- 
ing the damped system will have the identical form to the 
adiabatic equations, with the modification that the sound 
speed squared c^ is replaced by c^ . The damped form of the 
dispersion relation for an isothermal atmosphere is therefore 
taken as 



uj uJcayK^^ + (uj — a /C|| j X 
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= 0, (8) 



Pi 
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where a^ is the vertical component of the Alfven velocity and 
Oy is the component perpendicular to the plane containing 
k and g. N is defined by iV^ = g/H - g^/c^, <2jc = c/2H, 
and fch is the horizontal component of the wavevector. 

Although ad hoc, the above dispersion relation will 
prove useful in Section 3.1 in understanding the connectiv- 
ity of gravity waves to field guided acoustic (slow) waves or 
Alfven waves in the upper atmosphere where a ^ c. (See 
Newington & Cally (2010) for the low-/? asymptotic solu- 
tions of the 3D MHD dispersion relation.) As in paper I, 
we corroborate the dispersion diagrams with numerical so- 
lutions of the linearized wave equations, the derivation of 
which does not assume an isothermal atmosphere. This is 
described in the following section. 

2.2.4 Numerical solution of linearised, damped MHD 
equations 

The linearised, MHD equations incorporating Newton cool- 
ing are obtained from using the (non-isothermal) expression 
for the pressure perturbation (2), in the momentum equa- 
tion, and eliminating the density by means of the continuity 

-'- Note that this is equivalent to redefining the ratio of specific 
heats, as in Bunte & Bogdan (1994). 
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equation. Expressed in terms of the cartesian components of 
the Lagrangian displacement vector, ^ = (^, rj, Q the equa- 
tions are as follows: 



uj ^ — kx {igCm +c 112 {kx£, — iC')) 
a [k^ cos 9 cos sin 9 (^ + k^ cos (j) 
k^ cos 9 ^ ~ k^ sin 9 sin (f>£^ + ik., 
I kx cos 6 sin 9 sin (prj' — cos 9 cos sin ^ (^" + cos 9 ^"] = , 



+ a [k^ cos 9 cos sin 9 (^ + k^ cos (^ sin 6* sin (j) rj 
— k^cos OS^ — k^sin 9 sin ^ + ifca; sin 61 sin (/>(" 



(9) 



- -^M2 (ifc^ C + C') 

+ g ((1 - /ii) C' + ifc. - -|C (1 - Ml + ^m'i) 

+ (oj — kj^a cos sin 6 — Miff') C + i'^a; C Mac ' + /J2C '("' 

+ IKi C ^2? + IKi: C 5^2 + C C M2 + C ^2C 

+ a sin 6' [fc^, cos (j> cos ^ ^ — cos 61 (sin r;" + cos (f> ^") 

+ sin^ (ifea; sincjf) (sin(^^' — cos 77') + (^")] = , (10) 



(cj — fc^a cos (f> sin S) ?) 
+ a [sin6' (fcj; cos<^sin^sin(;/) (fcj,^ — i(^') 
+ cos 9 (2i kx cos 4>ri' ~ sin (ifca; ^' + ("")) ) + cos ^ j)"] = 



(11) 



The Alfven speed is denoted by a. 

The purely horizontal case 9 = 90° is not examined 
here as it is singular in nature, with the governing equations 
reduced in order resulting in 'critical levels' producing res- 
onant absorption at the Alfven and cusp resonances (Cally 
1984). This 'absorption' may in fact be interpreted as a mode 
conversion, where the absorption coefficient is continuous 
with 9 as it reaches 90° (as found by Cally & Hansen (2011) 
in the case of the Alfven resonance), so physically there is 
nothing particularly special about horizontal field despite its 
mathematical peculiarity. The singularities in the solutions 
are regularised in practice due to solar atmospheric waves 
not being strictly monochromatic. 

As in paper I, a driven wave scenario is envisaged, where 
monochromatic upward propagating (in terms of group ve- 
locity) gravity or acoustic waves are assumed to be excited 
at the bottom of our region of interest. In the interest of 
isolating the effects of mode conversion of acoustic and grav- 
ity waves during their propagation through the atmosphere, 
the condition that no slow magneotacoustic waves or Alfven 
waves entered the photosphere from below is imposed. Con- 
versely, no waves of any sort are allowed to enter from the 
top. 

The choice of a boundary value problem means that we 
considered the spatial damping of these waves (see Souffrin 
(1972)). Temporal invariance of the atmospheric model con- 
sidered means that the wave frequency lu is constant; this is 
determined by the driving force and is real. Horizontal in- 
variance ensures constancy of the horizontal component of 
the wave number kx during the wave's propagation. Damp- 
ing results in the 2 component of the wavenumber kz being 
complex. The magnitude of the imaginary component of kz 
increases with the strength of the damping. 

Solution of the equations requires an arbitrary normal- 
isation condition be applied. In paper I this was chosen to 



be normalisation of the vertical velocity, u; to 1 km^^ at 
the base of the photosphere (z — 0). For most of the re- 
sults presented in section 3, we chose to normalise w to 1 
km~^ at the base of the underdamped region (0.2 Mm), but 
we also ran simulations where the wave was normalised at 
the base of the photosphere (0 Mm) to observe the effect of 
overdamping on the transmitted fluxes. 

The motion is assumed adiabatic (r — > 00) at heights 
above 1.6 Mm, where the radiative decay time becomes very 
long, and below the height of application of the normalisa- 
tion condition. This simplification allows us to apply the 
same boundary conditions as in paper I, without modifica- 
tion. The reader is referred to that paper for details. 

The acoustic and magnetic fluxes were calculated in the 
same manner as in paper I. 



3 RESULTS AND DISCUSSION 

3.1 Dispersion diagrams 

We present diagnostic dispersion diagrams for four different 
orientations of the magnetic field. Figure 4 presents the adi- 
abatic results. This is the same diagram as figure 4 in paper 
I, except that in this paper the z axis has been truncated at 
the lower boundary of 0.2 Mm, and an essentially evanes- 
cent branch is included (yellow). Figure 5 displays results 
obtained for damped wave propagation, using the isother- 
mal dispersion relation, (8). Recall that the lower branch 
of the dispersion curve corresponds to the up-going gravity 
wave. The asymptotic solutions for the field guided acoustic 
wave and the Alfven wave are indicated in the region above 
the equipartition level, by the dashed blue and dotted red 
curves, respectively. 

In both figures, the curves are colour coded according 
to the damping per wavelength (Im {k}/ Re {fe}), where the 
more yellow the curve, the larger the value, and the heavier 
the damping. In the adiabatic results, (fig. 4), the yellow 
curves represent a predominantly evanescent mode that is 
absent from the figure in paper I. We draw attention to it so 
that its counterpart can be identified in the damped figure, 
but we are not concerned with this mode in this paper. The 
colour coding in figure 5 shows that the up-going gravity 
wave is most heavily damped lower in the atmosphere, which 
is to be expected from the linear form of r. 

The dispersion diagrams imply that the behaviours 
found in the adiabatic case, concerning the wave propaga- 
tion behaviour with the field orientation and the character 
of the dominant mode conversion, are largely unaffected by 
the damping. Comparison of figures 4 and 5 show that the 
wave paths in z — kz space are preserved when damping is 
included. This implies that the conclusions drawn in paper 
I about gravity wave behaviour with various field orienta- 
tions are also valid when the gravity wave experiences radia- 
tive damping. The connectivities to the asymptotic solutions 
(which are indicative of the dominant mode conversion) are 
preserved despite the introduction of damping. Numerical 
integration is used to confirm the dominant mode conver- 
sions in the following subsection. 

The paths in z-kz space are largely unchanged when 
constant values of r (ranging from 200s to Iks) are used 
instead of the linear form (5) . The strength of the damping 
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Figure 4. Dispersion diagrams for adiabatic propagation of waves of frequency 1 niHz and horizontal wavenumber kx=2 Mm~^, subject 
to a lOG magnetic field, included for the purpose of comparison to the radiatively damped results. In this dispersion diagram and those 
that follow, the vertical lines indicate the value of the ratio a?" ji? at various heights: solid line - a?" jc? = 1; dashed line - o? /c^ =0.1; 
dotted line - a?" jc? = 0.01; dot-dashed line — a? jc? = 10. The curve is coloured according to the damping per wavelength - the larger 
the magnitude, the greater the damping, and the more yellow the curve. The asymptotic solutions (detailed in paper I) corresponding 
to the field-guided acoustic waves and Alfven waves are shown as blue dashed, and red dotted curves, respectively. The bottom branch 
corresponds to the up-going gravity wave. The four panels differ in the field orientations - Top left panel: vertical field. The gravity 
wave is reflected as a slow wave; top right panel: 2D inclined field where the field lines are closely aligned with the direction of wave 
propagation. The gravity wave penetrates the equipartition level and undergoes mode conversion to a field-aligned acoustic wave; Bottom 
left panel: 2D inclined, anti-aligned field. When the attack angle (the angle between the direction of wave propagation and the magnetic 
field) is too large, the gravity wave reflects as a slow wave; Bottom right panel: 3D inclined field. Here the gravity wave connects to the 
Alfven wave solution. This figure presents the same results as figure 4 of paper I, but with the z axis starting at 0.2 Mm. Note also that 
the plots shown here include a predominantly evanescent mode (yellow) that is absent in the original figure. 



along the paths varies with each form of r. Of course, the net 
flux emerging at the top depends on this damping along the 
whole path, and is addressed in the next subsection. The 
results presented there also justify the crude assumptions 
made in our dispersion analysis ex post facto. 



3.2 Transmitted fluxes 

In all cases the radiative damping results in a reduction 
in magnitude of the measured acoustic and magnetic fluxes 



with respect to the adiabatic results. Performing the integra- 
tions with both linear and constant forms of r revealed that 
the magnitude of the flux reduction in the damped simula- 
tions is very sensitive to the form of the radiative damping 
used. Our investigations allow us to make qualitative state- 
ments regarding the effects of magnetic field strength and 
orientation, the effect of frequency, and the impact of the 
depth of the application of the normalisation condition, and 
we illustrate them with the figures obtained using (5). How- 
ever, bear in mind that the actual values of the measured 
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Figure 5. The same gravity wave propagation scenarios as in figure 4, except that in this case the propagation is subject to radiative 
damping. The colour of the curves indicates the damping per wavelength. Damping is greatest for the gravity wave in the photosphere, 
with the energy losses decreasing as the height increases. (Recall that the yellow curves with k^ tending to zero high up are the 
predominantly evanescent modes present in the adiabatic case.) Comparison with the adiabatic results (figure 4) illustrates that the 
mode conversion pathways found from the adiabatic simulations, are preserved in the presence of damping. 



fluxes in the figures are contingent on the use of the Mihalas 
& Toomre (1982) form of the radiative damping time. 



batic and damped cases. This was observed in simulations 
with constant r too. 



3.2.1 Effect of magnetic field orientation 

Figure 6, and the left hand panel of Figure 7 show the be- 
haviours of the magnetic and acoustic fiuxes as functions of 
the inclination 6 and azimuthal angle 4> for a gravity wave 
with frequency 1 mHz and horizontal wavenumber 2 Mm^^, 
in an atmosphere with an imposed magnetic field of lOG. 
Figure 6 shows the behaviour at three different azmiuthal 
angles. As expected from the dispersion diagrams, apart 
from the obvious reduction in flux magnitude, the behaviour 
of the fluxes with magnetic held orientation in the damped 
simulations is similar to that observed in the adiabatic case. 
Notably, the angles at which the ramp effect turns on and 
the mode conversion takes place are the same in the adia- 



3.2.2 Effect of magnetic field strength 

Figure 8 and the righthand panel of Figure 7 show the results 
calculated with a magnetic held strength of lOOG. Compar- 
ison with the results generated with a weaker fleld strength 
of lOG (flg. 6), reveals that the waves experience less flux 
attenuation in the high field strength cases. Although mag- 
netic pressure and tension play a larger role in driving waves 
at higher field strengths, this result is not obvious a priori. 



3.2.3 Effect of frequency 

Figure 9 shows the results calculated for two different fre- 
quencies. The top panel is for waves with frequency 0.7 
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Figure 7. Acoustic (full) and magnetic (dastied) wave-energy fluxes (kW m ^) as functions of magnetic field orientation with inclination 
85° . Results of damped and adiabatic simulations are the black and grey curves, respectively. 



mhz, and the lower panel is for 2.1 mHz. This shows that 
the low frequency waves experience greater flux attenuation 
than waves of higher frequency. This is to be expected since 
shorter period waves naturally lose less energy radiatively 
at fixed r. 



3.2.^ Effect of the depth of application of normalisation 
condition 

Figure 10 shows the fluxes as a function of magnetic field 
inclination generated with the normalisation condition ap- 
plied at the base of the photosphere at 2 = Mm. In these 
simulations, the gravity wave is effectively extinguished, as 
was predicted by Mihalas & Toomre (1982). 



4 CONCLUSIONS 

The main conclusions we draw from the results of our inves- 
tigation are: 

(i) The primary effect of radiative damping in our model 
is a decrease in the measured wave energy fluxes obtained, 
with respect to the adiabatic results. The extent of the 
flux attenuation is highly sensitive to the radiative damp- 
ing times. 

(ii) Damping does not have a great effect on the mode 
conversion pathways in the z — kz plane. Speciflcally, direct 
coupling to Alfvenic disturbances higher up in the chromo- 
sphere seems to survive (in diminished form) the introduc- 
tion of radiative losses. 

(iii) Short damping times very effectively extinguish grav- 
ity waves. This is most likely to be a problem in the low 
photosphere, where radiative damping times are expected 
to be far shorter than wave periods, rendering the waves 
overdamped. Reconciling observations of propagating pho- 
tospheric gravity waves (Straus et al 2008; Stodilka 2008) 
with this seems an impossibility unless they are being gener- 
ated there in situ? Komm et al. (1991) suggest that gravity 

^ In Earth and planetary atmospheres, it is known that primary 
gravity waves can produce secondary gravity waves by nonlinear 



waves are generated by decay of granular motions in the mid 
photosphere. Stodilka (2008) found that overshooting con- 
vection starts at around 100 km and apparently detected 
gravity waves above this height. Gravity wave excitation by 
penetrative convection (albeit in the context of the convec- 
tion zone/radiative core interface) has been modelled numer- 
ically (Dintrans et al. 2003, 2005) and found to be effective. 
(iv) Even accepting in situ gravity wave generation, it is 
not clear that Straus et al (2008) see propagating gravity 
waves as low as 70 km, where they should be substantially 
overdamped, notwithstanding the rather arbitrary form of 
our radiative decay time t{z). Indeed, their own numerical 
simulations using the (nonmagnetic) radiation hydrodynam- 
ics CO^BOLD code, which has a much more sophisticated 
treatment of radiation than we have instituted, is dominated 
by convection at ~ 70 km, with a negligible gravity wave flux 
there (Straus, private communication). Gravity waves are 
progressively generated as height increases, reaching their 
full strength at around 300 km, with the gravity wave being 
'free' only above this. This justifies our decision to model 
gravity waves in 2; > 200 km; the region below 200-300 km 
is where the waves are driven. 
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Figure 9. Acoustic (full) and magnetic (dashed) damped (black) 
and adiabatic (grey) wave-energy fluxes (kW ia~-^) as functions 
of magnetic field inclination 8 with (p = 30° for B=10 G, kx = 2 
Mm~^ for waves of two different frequencies. Top panel: 0.7 mHz; 
Lower panel: 2.1 mHz. Note the different 6 scales on the two 
graphs. 
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Figure 10. As for 6, but with the normalisation condition now 
applied deeper at z=0. Very little flux is measured at the top 
when the gravity wave is damped. With the chosen form of t{z), 
the short damping times in the low photosphere very effectively 
mitigate gravity wave progression. 



